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The existence of discrete breathers (DBs), or intrinsic localized modes (localized periodic oscillations of 
transzigzag) is shown. In the localization region periodic contraction-extension of valence C — C bonds oc- 
curs which is accompanied by decrease-increase of valence angles. It is shown that the breathers present in 
thermalized chain and their contribution dependent on temperature has been revealed. 



I. INTRODUCTION 

Localized excitations in nonlinear systems (solitons, po- 
larons, breathers) became the subject of growing interest dur- 
ing last decades. Discrete breathers, or intrinsic localized 
modes correspond to shortwave length periodic vibrations 
which can not be revealed in continuum model of the chain. 
Intensive study of the breathers has an origin in pioneer work 
of Sievers and Takeno [1]. Their existence is proved by cor- 
responding theorem [2, 3] and numerous numerical investiga- 
tions [4] . Now their role in mechanisms of the energy relax- 
ation in molecular systems becomes rather clarified [5, 6, 7]. 

However all these results relate to simple one-dimensional 
systems. Meantime it should be waited that similar elemen- 
tary excitations can exist in more realistic discrete nonlinear 
models [8]. The necessary conditions for their existence are 
discreteness of the system leading to boundedness of phonon 
frequency spectrum and amplitude dependence of eigenfre- 
quencies caused by nonlinearity. These conditions are satis- 
fied for polymeric molecules. For example, anharmonic po- 
tential of C — H valence bond in carbon-hydrogen chains leads 
to existence of high-frequency localized vibrations of such 
bonds [8]. Naturally, this type of breathers has to be present 
in polyethylene (PE) chain but corresponding frequencies are 
very high (~ 3100 cm -1 ). We will discuss the breathers 
caused by the backbone vibrations of PE chains which cor- 
respond to much more less frequencies because they arise 
due to coordinated changes of C — C valence bonds and CCC- 
valence angles. 



n. MODEL OF MACROMOLECULAR CHAIN IN PE 
CRYSTAL 

Nonlinear dynamics of planar PE macromolecule has been 
considered in [9, 10] where a detail consideration of the tran- 
szigzag model is presented. This is a reason for only short its 
description which is given below 

PE crystal consists of parallely packed zigzag chains [PE 
macromolecules (CH^ooL The structure of the crystal is 
schematically plotted at Fig. 1. Every PE macromolecule 
in the crystal has a transzigzag conformation. It means that 



backbone of the chain has a planar zigzag structure with equi- 
librium length po = 1.53 A of valence bond H2C-CH2 
and value 9 = 113° of valence angle CH 2 -CH 2 -CH 2 . All 
macromolecules are situated in parallel planes. The crys- 
talline structure is characterized by angle 8 and three periods: 
a = 4.51 A, b = 7.031 A, c = 2p sin(0 o /2) = 2.552 A . 




0.5 b 



-0.5b 



FIG. 1: Schematic representation of the crystalline PE. The central 
trans-zigzag backbone (curve 0) and the six neighbor chains (curves 
1,2,. ..,6) are shown. The local coordinates for the central chain are 
presented. 

Let us consider the dynamics of PE macromolecule and 
will take into account its interaction with six neighbor chains 
which being considered as immobile (approximation of im- 
mobile neighbor chains). Because the motion of hydrogen 
atoms with respect to backbone chain is not significant when 
studying low energy nonlinear dynamics, we will consider ev- 
ery CH2-group as a united particle. 

Let us introduce the local coordinate system presented at 
Fig. 1. Then the Hamiltonian of the chain can be written as 
follows 
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where the first term describes the kinetic energy of nth unit, 
the second the deformation energy of the nth valence bond, 
the third the deformation energy of the nth valence angle, the 
fourth the deformation energy of nth torsional angle, and the 
last term the energy of interaction of the nth unit with the six 
neighboring chains (substrate potential), M = 14m p is the 
mass of the united atom (m p is the proton mass). 
The length of nth valence bond is 

/ 2 ,2 i 2 \l/2 
Pn = K,l + a n,2+ V3) > 

where a n>1 = u n + u n +i - l x , a n.2 = v n + v n+ i, a n , 3 = 
w n+1 - w n + l z [l x = p cos(6» /2) and l z = p a sin(0 o /2) 
- transversal and longitudinal step of the zigzag chain]. The 
cosine of the n-th valence angle is 

COs(0 n ) = (a n -l t ia n ,l + a n -l,20.n,2 ~ On- 1,30x1,3) I Pn-lPn, 

and the cosine of the ?ith torsional angle is 

C0S(5„) = (— &„,l&n+l,l - &n,2&n+l,2 + K,sK+l,3) /finfin+1, 

where 

b n ,l = an-l,20n,3 + a n,2 a n-l,3j 

b n ,2 — an-l,lOn,3 + a n,\0-n-l,Zi 

bn,3 = Cn-l,2&n,l — 0'n,20"n-\,li 

Pn = K 1+ bl 2 + bl 3 ) 1/2 . 

Potentials of valence bond, valence angle, and torsion angle 
are accepted in the form 

V(p n ) = D {1 -cxp[-a(p- po)}} 2 , 

U(0n) = ^(cOsfl-COS^o) 2 , 

W(S n ) = C\ + C 2 cos S n + C 3 cos 35 n , 

where the parameters Do — 334.72 kJ/mol, a = 1.91 A~\ 
7 = 130.122 kJ/mol, d = 8.37kJ/mol, C* 2 = 1.675 kJ/mol, 
C3 = 6.695 kJ/mol. The substrate potential 

Z(u, u, w) = e w sin 2 (7rw/Z z ) 
+\K U [l + £ u s\vl 2 (ttw/I z )] {u - \l x [1 - cos(irw/l z )}} 2 
+\K V [1 + e v sm 2 {irw/l z )] v 2 , 

where the parameters e u = 0.0674265 kJ/mol, e v = 
0.0418353 kJ/mol, e w = 0.1490124 kJ/mol, K u = 2.169513 
kJ/Amol 2 , K v = 13.683865 kJ/Amol 2 . Detail substantiation 
of chosen potentials is presented in [10]. 



III. LINEAR OSCILLATIONS OF TRANSZIGZAG 

Small-amplitude vibrations of isolated transzigzag were 
considered in [9, 11, 12] and vibrations with account of in- 
teraction with immobile surrounding chains - in [10]. 

Small-amplitude vibrations can be divided on planar (in the 
zigzag plane) and transversal ones. In turn, the plane motions 
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FIG. 2: Dispersion curves u> = 0J a (q), IjJ = Wt(g), w = co (q) 
(curves 1, 2, 3) for trans-zigzag interacting with immobile surround- 
ing chains (a). Density of energy distribution p on frequencies u> for 
thermal vibrations with temperature T — IK (b), T — 100 K (c), 
T — 200 K (d), and T — 300 K (e). Gray color corresponds to 
frequencies region in which the discrete breathers occur. 



are divided on low-frequency acoustic and high-frequency 
optic vibrations. Accordingly, one can separate three dis- 
persion curves corresponding to 1) plane acoustic phonons 
u> = u> a (q)> 2) plane optic phonons uj = uj (q), 3) transversal 
phonons to = u>t(q). These curves with account of interchain 
interaction are presented at Fig. 2 (a). 

For isolated chain (substrate potential Z(u, v, w) = 0) the 
acoustic phonons have the frequency spectrum < uj a < 
228.75 cm -1 , torsional (transversal) phonons - spectrum < 
u>t < 243.85 cm -1 , and optic phonons - spectrum 838.47 
cm -1 < uj < 1168.11 cm -1 . Interaction with surrounding 
chains leads only to a shift to right of low boundary of acoustic 
and transversal phonons: 19.05 cm^ 1 < uj a < 229.54 cm -1 ; 
52.43 cm" 1 < uj t < 249.43 cm" 1 ; 838.69 cm" 1 < uj d < 
1168.27 cm- 1 . 



IV. LOCALIZED NONLINEAR VIBRATIONS OF PE 
CHAIN 

The equations of motion corresponding to Hamiltonian (1) 
have the form 
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Complexity of this system does not allow to find its analytic 
solution. Thefore we used a numerical procedure. 

The finite chain consisting of N — 200 CH2 groups was 
considered. The viscous friction providing adsorption of 
phonons was introdused on the boundaries of the chain. The 
system of equations (2) with n = 1, 2, N was integrated 
numerically with breather-like initial condition. If a discrete 
breather can exist in the chain it can manifest itself in irra- 
diation of superfluous (non-breather part of initial excitation) 
phonons. 

Numerical modelling has shown that only one type of local- 
ized periodic vibrations exists caused by tension-compression 
of valence C — C bonds with coordinated change of valence 
angles in the plane of transzigzag - see Fig. 3 (b) (c). The 
vibration occurs in the plane of trans-zigzag with nodes dis- 
placements perpendicular to main backbone axis [Fig. 3 (a)]. 

These vibrations are stable excitations, which are charac- 
terized by frequency u>, energy E and dimensional width 
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where the point n c = J2 n n Pn determines the position of the 
vibrations center and sequence p n = *}2 n E n /E - the dis- 
tribution density of energy along the chain. Essential non- 
linearity of these vibrations is manifested in decrease of its 
frequency with amplitude growth. So, the revealed excitation 
is actually discreet breather. 

Dependence of energy E and width of the breather L on 
its frequency u> is presented at Fig. 4. The frequency spec- 
trum of the breather has situated near low boundary of optic 
phonons. It is a reason why weak intermolecular interaction 
with surrounding chains do not practically effect on such ex- 
citations. If in isolated chain its frequency uji, = 817 cm" 



< 

in the chain surrounding by immo- 

u b = 820 cm" 1 < u> < w o (0) = 839 



ui < w o (0) = 838 cm" 1 
bile neighbor chains 

cm -1 . When frequency decreases the energy of breather in- 
creases monotonically and its width decreases monotonically. 

The frequency of breather has to be separated from the fre- 
quencies of small amplitude linear vibrations. One could ex- 
pect also the existence of low-frequencies breathers with fre- 
quencies near upper boundary of transversal phonons ui > 
maxw t (q) and high-frequency breathers with frequencies 
near upper boundary of optical phonons ui > u) (tt). But nu- 
merical analysis has shown that such localized excitations are 
absent. 
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FIG. 3: Localized planar periodic vibrations of trans-zigzag . Vibra- 
tions are schematically shown, the thickness of line corresponds to 
amplitude (a). The magnitudes of valence bonds p n (b) and angles 6 n 
(c) are presented for ten different instants. The frequency of breather 
u = 820.5 cm -1 , energy E = 26.4 kj/mol, width L = 4.28 



V. THERMAL VIBRATIONS OF TRANS-ZIGZAG AS 
ORIGIN OF DISCRETE BREATHERS 

Let us consider the thermal vibrations of trans-zigzag. With 
this goal we analyze the finite chain consisting of N segments. 
Their No segments near boundary (from both sides) are situ- 
ated into heat bath with temperature T. The dynamics of the 
system is described by the systems of Langevin equations 
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+ Cn - T n Mu n , 

+ Vn- r n Mu„, 
+ Cn - r„Mu„, 

n = 1,2, ...,N , 



(3) 



where the Hamiltonian of the system H is given by Eq. (1), 
£„, rjn, and £„ are random normally distributed forces describ- 
ing the interaction of nth molecule with a thermal bath, the 
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FIG. 4: Dependence of energy E (a), width L (b) of the breather 
upon frequency u in isolated chain (curves 1, 3) and in the chain 
with substrate potential (curves 2, 4). 

coefficient of friction r„ = for Nq < n < N — Nq and 
r„ = T for n < N and N - N < n < N. Coefficient of 
friction T = 1 /t r , where t r - the relaxation time of the veloc- 
ity of the molecule. The random forces £ n , r\ n , and £„ have 
the correlation functions 

(Zn(tl)U{t2)} = (Vn(tl)Vm(h)) = (Cn(*l)Cm(*2)) = 

2MTk B TS nm S{t 1 -t 2 ), 

(€n(*l)»7m(t2)> = <^n(tl)Cm(*2)> = <^(*l)Cm(ta)> = 0, 
1 < n, m < N , N - N < n,m < N, 

where ks is Boltzmann's constant and T is the temperature of 
heat bath. 

The system (3) was integrated numerically by the standard 
forth-order Runge-Kutta method with a constant step of inte- 
gration At. Numerically, the delta function was represented 
as 5{t) = for \t\ > At/2 and S(t) = 1/At for \t\ < At/2, 
i.e. the step of numerical integration corresponded to the cor- 
relation time of the random force. In order to use the Langevin 
equation, it is necessary that At <C t r . Therefore we chose 
At = 0.001 ps and the relaxation time t r = 0.1 ps. 

Let us consider a frequency distribution of kinetic energy 
of thermal vibrations. For this goal the system (3) was inte- 
grated numerically for N = 500, TVo = N/2. While choosing 
the initial conditions as corresponding to ground state of the 



chain the system was integrated during t = I0t r to bring it 
in the thermal equilibrium. After that we calculated the den- 
sity of molecules kinetic energy distribution on frequencies 
p(u>). To increase an accuracy, the density of distribution was 
calculated using 1000 independent realizations of the chain 
thermalization. The profile of the density of distribution for 
different values of temperature is presented at Fig. 2 (it is 
accepted that J p(Ld)duj = 3). 

For temperature T = 1 K the density of distribution co- 
incides practically with corresponding density for linearized 
system, so anharmonicity is here not essential. All vibra- 
tions are linear and only phonons are thermalized. With fur- 
ther increasing of temperature the amplitude of thermal vibra- 
tions also increases so their anharmonism is manifested. For 
T = 100 K we see a shift of density behind the low boundary 
of the spectrum of optic phonons which becomes more pro- 
nounced with further increase of the temperature. High fre- 
quency vibrations in this region can be identified as breathers 
because they arise in the frequency region [uib, w o (0)], corre- 
sponding to this type of excitations. The part of energy corre- 
sponding to breathers may be found as 

/•o)o(O) 

Pb = p(u)duj. 

The contribution of the breathers in thermal energy increases 
with growth of temperature (for T = 1 K it is pj, = 0.002, 
for T = 100 K — Pb = 0.106), has a maximal value pb 
0.1 IT, for T = 200 K and then decreases (for T = 300 K the 
breathers contribution corresponds to p = 0.083). 

Let us isolate the breathers from thermal vibrations. With 
this goal we consider the chain consisting N = 500 segments 
with boundary cites (AVj = 50) embedded to the heat bath 
with temperature T. After thermalization, we put the temper- 
ature of thermal bath T = and consider the irradiation of 
heat energy from internal region (A*^ < n < N — No). The 
relaxation process for T = 200 K is shown at Fig. 5. We see a 
formation of several mobile localized excitations. Their detail 
analysis leads to conclusion that they are discrete breathers 
with frequencies oj ~ uj,. So, one can observe the presence of 
breathers in thermal vibrations. 

Let us consider the interaction of discrete breather with 
thermal phonons. The stationary discrete breather with fre- 
quency oj = 820.5 cm -1 was situated in the center of finite 
chain (N = 200), their edges (No = 10) being situated in 
thermal bath (T = 10 K). As we can see from Fig. 6 the 
breaking of breather is observed just as the center of the chain 
is thermalized (the energy loss is 50% for 20 ps). 

Probability of thermally activated formation of discrete 
breathers in the chain grows with increasing of the tempera- 
ture. Therefore their concentration has to increase when tem- 
perature grows. However in a thermalized chain the breather 
has a finite time of life, decreasing with growth of the temper- 
ature. It is reason for non-monotones dependence of the con- 
centration of breathers pj, upon temperature T — it increases 
when T < 200 K and decreases for T > 200 K with maxi- 
mal magnitude for T = 200 K. Numerical study shows that 
the breathers may be better separated from thermal vibration 
namely when T = 200 K. 
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FIG. 5: Formation of discrete breathers from thermal vibrations of 
zigzag chain (N = 500, T — 200 K). The absorbing ends are 
considered (No = 50). Temporal dependence of energy distribution 
E n in the chain is presented. 




FIG. 6: Breaking of the breather in thermalized chain (N — 200, 
iV = 10, T = 10 K, frequency of the breather lo = 820.5 cm" 1 ). 
Temporal dependence of current local magnitudes of temperature (ki- 
netic energy of chain segments) T n is presented. 



One can see from Fig. 2 that revealed breathers is a 
unique type of stable localized periodic excitations in ther- 
malized chain for given parameters of the crystal. Besides the 
breathers only vibrations with frequencies possessing to spec- 
trum of linear oscillations can be thermalized. It confirms the 
conclusions with respect to existence of only one stable type 
of discrete breathers corresponding to localized oscillations 
of valence C — C bonds. The breathers present in thermalized 
chain even for sufficiently small temperatures. 



VI. CONCLUSION 



Stable localized nonlinear vibrations which are discrete 
breathers can exist in polymer crystal. In PE macromolecule 
they are planar vibration of transzigzag with periodic defor- 
mation of valence bonds C — C and valence angles CCC. The 
breathers present in thermalized chain and their contribution 
in heat capacity may be essential. 
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